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We study the electromagnetic properties of a hot medium at temper- 
atures below electron mass. It was observed earlier that the first order 
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^\ , Abstract 
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^— ^ ' hot loop corrections do not affect the electromagnetic properties of hot 

\ jr^ , media due to the absence of self-interaction of photons. However, the 

f~^ ' second order contributions are found to modify these properties due to 

the overlapping hot and cold loops. It has also been noticed that the hot 
loops have to be integrated before the cold ones to establish the order by 
order cancelation of singularities to ensure renormalizability of the theory. 
CL(, It is also explicitly shown that the perturbative series of QED is still a 

D . convergent series in a covariant formalism at low temperatures. 



1 INTRODUCTION 

In quantum field theory, the thermal background effects are incorporated through 
the radiative corrections. The renormalization of gauge theories at finite tem- 
perature and density (FTD) requires renormalization of gauge parameters of the 
corresponding theory. Propagation of particles and the electromagnetic prop- 
erties of media are also known to change in FTD framework and this change 
could be determined through renormalization techniques. The massive parti- 
cles thereby get enhancement in their masses at one-loop [1-3], two-loop[4] and 
presumably at all loop levels. However, the gauge bosons acquire a first order dy- 
namically generated mass due to plasma screening efFect[5,6] in gauge theories. 
It helps to determine the changes in electromagnetic properties [7] of a medium. 
This effect has already been studied in detail at the one-loop level and showed 
it explicitly [5] that the electric permittivity and the magnetic susceptibility of a 
medium are not modified by incorporating the first order thermal background 
effects. At low temperatures, i.e. T « We (where rrie is the electron mass). 



the vacuum polarization tensor in order a does not acquire any hot corrections 
from photons in the heat bath because of the absence of self-interaction of pho- 
tons in QED. However, the higher order corrections to the vacuum polarization 
tensor of photons are nonzero in the same background. The electromagnetic 
charge in QED modifies in such a situation and thus leads to the modifications 
in the electromagnetic properties of a medium. These statistical contributions 
are calculated either in Euclidean or Minkowski space by using imaginary or 
real-time formalism, respectively. In Euclidean space the covariance breaks and 
time is included as an imaginary parameter. On the other hand, in the real-time 
formalism, an analytical continuation of energies along with Wick's rotation re- 
stores covariance in Minkowski space at the expense of Lorentz invariance. It 
has been noticed earlier that the breaking of Lorentz invariance can lead to 
the non-commutative nature of gauge theories [8]. This covariance is incorpo- 
rated through 4-component velocity of background heat bath which is given as 
u^ = (1, 0, 0, 0). The Lorentz invariance breaking hot term and conserving cold 
terms remain separate at the one loop level since the propagators comprise of 
temperature dependent (hot) terms added to temperature independent (cold) 
terms. This is how the simplification in calculations is ensured by segregation 
of hot and cold terms in this formalism. At the two- loop level an overlap of hot 
and cold terms occurs and the specific order of integration is the only way to 
get correct finite results. At the higher loop level, this combination (overlap) of 
hot and cold terms in the propagators makes the situation somewhat compli- 
cated and one needs to use some specific techniques to get rid of singularities. 
At low temperatures, 6{0) type pinch singularities appear in Minkowski space. 
This problem has been resolved in thermoficld dynamics [9] by doubling the 
degrees of freedom. Therefore, the particle propagators become 2x2 matrices 
whose 1-1 elements correspond to the usual thermal propagators [1], i.e., in case 
of fermions 

S0{p) = {^-m) {-, '-^—- - 27tS{p^ - m^)nF{Ep) \ , (1) 



p^ _ j^.i ^ jg 



where 



--(^^) - -^^i^i ' (2) 



is the Fermi-Dirac distribution function with (3~^. The boson propagator is 
taken as 

^r (P) = ^Vr-- - 2^Kk')nB (k)] (3) 

with 

""BiEk) - ^^(J)_^ , (4) 

is the Bose-Einstein distribution function in the usual notation. The rcnormal- 
ization of QED in this scheme [1-5] has already been checked in detail at the one 
loop level for all temperatures and chemical potentials. At the higher- loop level. 



the loop integrals involve a combination of cold and hot terms which appear due 
to the overlapping propagator terms in the matrix element. These higher or- 
der integral calculations are undoubtedly much more involved and some special 
techniques are needed to be used for them up to the two loop level. Even higher 
loops may not be manageable analytically. The electron self-energy has already 
been calculated at the two-loop level [4] . The photon self-energy calculations are 
even more involved so we restrict ourselves to the low-temperatures to explicitly 
prove the renormalizability of QED at the two-loop level through the order by 
order cancelation of hot singularities. In the next section we give the detailed 
calculations of the vacuum polarization tensor. Section 3 is comprised of the 
calculations of the renormalization constants of QED up to the two-loop level. 
The electromagnetic properties of the hot medium are discussed in section 4. 
Discussions of these results is included in section 5. 



2 VACUUM POLARIZATION TENSOR IN QED 



The vacuum polarization tensor of photon at the two-loop level gives the second 
order hot corrections to the charge renormalization constant of QED at low 
temperature. This contribution basically comes from the selfmass and vertex 
type of electron loop corrections inside the vacuum polarization tensor. The 
diagrams in Fig.l give the main thermal contribution through the hot virtual 
photon to the vacuum polarization tensor up to the order o? at low temperature, 
i.e., T <C vTig. In our scheme of calculations, the self energy of photon in Fig. (la) 
is given by 



n;i.(p) = n^.(p,r = o)- 



while that in Fig. (lb) is 



n^.b) = n^.(p,r = o)- 






d*k 



dH 

(2^ 



dH 



(27r)4 J (27r) 



(5) 

■J^,Sf3{kyfpD''p'{l)Sf}{k+lyf,SfJ{k+l-pyfpSp{k-p). 

(6) 



The hot terms give the overlapping divergent terms and the removal of this 
divergence to establish the renormalization is done using different calculational 
techniques, which includes thcrmofield dynamics also. We have used an alter- 
native technique to get rid of these (5(0) type pinch singularities [10] through 
the identity. 



[Sik'f 



5{e 



1 



5\k\ 



(7) 



B + ie 2 

This corresponds to the evaluation of the 1-1 component of the matrix prop- 
agator only. Doubling the degrees of freedom is not required in this approach 
because the order by order cancelation of singularities can explicitly be shown 



by the use of Eq.(7). In this way one gets rid of the pinch singularities without 
summing over the unphysical ghost terms of hard-loops. However, in this tech- 
nique, the accuracy of results depend on the order of integration. One has to 
use the right order of integration to get the correct results. The renormalizabil- 
ity can also be established for the right order of integration only. Integration 
of thermal loop-integrals has to be done before the temperature-independent 
cold-loop integrals. The justification of this order could be found in the fact 
that the hot contribution in this technique correspond to the contribution of 
real background particles on mass-shell and incorporates thermal equilibrium. 
However the breaking of Lorentz invariance changes these conditions for the 
cold integrals. So the rcnormalization can only be proven with the right order 
of integration. To establish the renormalizability, it is sufficient to include the 
1-1 component of the matrix propagator only. We do not have to include the off- 
diagonal elements of the propagator matrices in our calculations. Moreover, the 
other components of the matrix propagators can be related to 1-1 components 
in thcrmoficld dynamics [9]. We restrict ourselves to low temperatures where 
hot fermion contribution in background is suppressed and the only thermal con- 
tributions come from the hot photons in the background heat bath. Therefore, 
the cold fermion and hot photon propagators are included. The calculations 
are simplified if the temperature dependent integrations are done before the 
temperature independent ones. The cold loops can later be integrated using 
the standard techniques of Feynman paramctrization and dimensional regular- 
ization, as is done in vacuum and are discussed in standard textbooks [11] . 
As an illustration of the importance of the correct order of integration, we just 
compare the results of one term. Let's consider the singular terms when the hot 
loop in Fig. (lb) is evaluated before the cold one, we simply get 

.9-H;:,(p,T) = ^-(^1--J (8) 

whereas, in the same term, the evaluation of the hot loop after the cold one 
gives 

9-K.iP^ T) . --,^- - ^ - ^C (3) (3|p| + -p„ j + ^. (9) 

The details of the calculations for Eqs.(8) and (9) are given in Appendix. This 
order of calculations has to be maintained everywhere. Same order of integra- 
tion was also maintained in the calculations of electron self-energy in Ref. [4]. 
This was the only way to show the order by order cancelation of singularities 
in the calculation of the second order electron mass even at high temperatures, 
i.e; T^ rrif.. The calculation of thermal corrections to the self-energy of pho- 
ton becomes very cumbersome at high temperatures. However, it would help 
to determine the corresponding changes in the electromagnetic properties of 



hot media. Although it might be very long calculations and difficult to do an- 
alytically without using any approximations, especially at temperatures near 
electron mass. The calculation of contribution of density effect will also be 
pretty complicated at the two-loop level. However the first order contributions 
have already been estimated for dense media in this scheme[3,5]. 

In a covariant formalism, the physically measurable couplings can be eval- 
uated through the contraction of the vacuum polarization tensor II^^, with the 
metric in Minkowski space 5"^" and u'^w [5]. This helps in evaluating the 
longitudinal and transverse components of vacuum polarization tensors in order 
to evaluate the thermal contribution to the electric permittivity and magnetic 
susceptibility. The gauge invariant finite contributions of the sum of all the 
diagrams of Fig.(l) are 



u^u^Tl,MT) = -'-^(l + ^,], (10) 



pI 



and 



5-n^.b,T) = ^. (11) 



Moreover, the longitudinal and the transverse components of the vacuum 
polarization tensor that can be calculated from Eqs.(lO) and (11) are 

n.(p,T) . ^^u^un,AP,T) ^ ^ (1 + ^) , (12) 

and 



IiT{p,T) = -hnUp^T) - g^-'UUp,T)] "'^' 



2 |p|2 i^ ^ 2m2 



respectively. These components of the vacuum polarization tensor can then 
be used to determine the electromagnetic properties of a medium with hot 
photons and will be discussed in the next section. It is worth-mentioning that 
Fig. (lb) gives major contributions to the finite terms in the above equations. 
However, the other two diagrams are needed to establish the cancelation of 
singularities even at low temperature. It is also interesting to note that the order 
of singularities is the same as the order of a in the matrix element. It is not an 
accidental occurrence of the matching order of singularity. It is due to the fact 
that the number of loops and the number of maximum delta functions in a term 
in the matrix element are the same and would correspond to the highest order 
of hot divergences. In the next section we have evaluated the renormalization 
constants of QED up to the two loop level at finite temperature. 



3 THE RENORMALIZATION CONSTANTS 

In this section we discuss the thermal contributions to all of the rcnornializa- 
tion constants of QED namely electron mass, electron charge and the electron 
wavefunction. Electron charge renormalization gives the thermal corrections to 
the coupling constant of QED, i.c; the fine structure constant. 



3.1 THE ELECTRON MASS 

The second order in a corrections to the electron self energy at low temperature 
have been calculated [4] to be equal to 



E0(p) = E^(p,T = O)- — 



where the singularity 



i[(i) + 6m)lA - 21] - 3(i. + Am)lA - ^1^.I'' - ^^(i^ - m) 

(14™ 

- = --7 -In 2~ ' (15) 

e rj TO^ 



is the usual ultrviolet divergence in the vacuum field theory in MS bar scheme 
of renormalization with 7 as the Euler- Mascheroni constant and 

^flBik), (16) 

/" 2tt^T\ 1+v 

— = TT- In , (17) 

E 3E^v 1-w' ^ ' 



I.p 27T^T^ r, ^ 1 

T 



2-,^2,„.[l-^--2-], (18) 



with V = — [2,12]. The divergences of the form Ia and - in Eq.(14) cancel on 
addition of the appropriate counter terms. The finite contribution then leads to 
the physical mass of electron: 



2 _ 2 



phys — '"■ l-"^! o 2 o 

'^ " \ 6m 6m 



(19) 



The electron mass and charge renormahzation can be obtained from the 
two loop self-energy for electrons and photons respectively. From Eq. (19) the 
change in electron mass due to temperature up to the order a^ relative to the 
cold electron mass m is 



8m _ T^ 
m rn? 



[a-K + A,o?) . (20) 



The renormalizability of the sclfmass of electron is proven not only through 
the order by order cancellation of singularities at both loop levels but also 
through the convergence of the pcrturbative series. It can also be seen in Eq.(20) 
that the second order term is much smaller than the first order term in the above 
equations. However, the situation looks a little different in charge renormahza- 
tion because of the absence of the first order thermal corrections to the photon 
self-energy at low temperatures. 

3.2 CHARGE RENORMALIZATION 

Using Eq. (13) from the last section and the standard method of the evaluation 
of the charge renormahzation constant of QED, the electron charge renormal- 
ization up to the order o? can be expressed as 

^3 = l + ^r^. (21) 

One can easily notice the absence of the first order temperature dependent 
term in Eq.(21). This is due to the absence of self interaction of photons in 
QED and has been discussed in detail in Ref. [5] . The second term in the above 
equation, the second-order contribution is the vacuum term is much smaller 
than 1 (vacuum part). The corresponding value of the QED coupling constant 
comes out to be 

afl = a(T = 0)(l + — ^). (22) 

It can be clearly seen in the above equations that the electron charge and hence 
the QED coupling constant goes smaller with the increased order of loops which 
assures the renormahzation of electron charge clearly. This change in the cou- 
pling constant leads to changes in the electromagnetic properties of a medium 
which we will discuss in the next section. 

3.3 WAVEFUNCTION RENORMALIZATION 

Using the same standard methods of QED and using the electron charge renor- 
mahzation constant, the wavcfunction renormahzation constant comes out to 
be 






3 + 1 






(23) 



It is noticeable again that the higher order terms are dominated by the lower 
order terms which assures the renormalizability of QED upto the two-loop level 
at finite temperature. In the next section we discuss the corresponding changes 
in the electromagnetic properties of a medium at low temperature. 

4 ELECTROMAGNETIC PROPERTIES IN A 
MEDIUM 

This is a well-known fact that the electromagnetic properties of media depend 
[7] on its statistical properties. However one-loop corrections does not change 
it much at low temperature whereas the second order corrections contribute to 
them. It can be seen using Ref. [5] that the electric permittivity of the medium 
at the two loop level modifies to 



eip,T) ^ 1 - \p\^IIl{p,T) ^ 1 



2ti2^2 



2a^T'p 



2m2 



(24) 



whereas the magnetic susceptibility is given by 



1 



^^ip,T) 



= 1- 



1 

9 



vl 



[nT(p,T)-f^ni(p,T)] = i- 



a'T' 



1 



pI 

2inn? 



(25) 

Equations (24) and (25) show that the second order radiative emission and 
absorption of hot photons from the heat bath leads to deviations from unity 
in the values of the dielectric constant and the magnetic permeability. This 
particular feature is not present at the one-loop level at low temperature because 
of the absence of FTD corrections to 11^ and Ht for cold fermions. There are 
two possible limits for 11^ and Ht as p^ — >0 and their physical interpretations 
depend on the order of taking the limits p^ — and |p| — > 0. In Eqs. (21) and 
(22), if po = IpI in the rest frame of the heat bath and then the limit |p| — > 
is taken, we get 



lim 0^(0, IpI, T) 

IpI^o 



0, 



(26) 



and 



'^pI 



t\in~' 



hm nT(0,|p|,T) = -— . 
IpI^o 2 



(27) 



On the other hand if we set pa = with |p| — > then we obtain 



coi^ hm ni(|p|,|p|,T)-0, (28) 

IpI — '0 

^1^ hm nT(|p|,|p|,T)-— ^. 
IpI^o 6 

This change in electromagnetic properties have been discussed in a httle more 
detail in the next section and the possible implications of these results are also 
mentioned thereafter. 

5 RESULTS AND DISCUSSIONS 



It is known from Refs.[5] that the low temperature effects on vacuum polar- 
ization in the photon background at the one loop level are zero because of the 
absence of self couplings of photons. At the higher loop level the vacuum polar- 
ization contribution is non-zero, even at low temperature due to the overlapping 
hot and cold terms. Most of these terms are finite or simple to regulate through 
the addition of all of the same order diagrams. However, the complications arise 
due to the pinch singularities in overlapping divergences. The pinch type sin- 
gularities are treated using the identity in Eq.(7). The renormalization of the 
theory can only be proved if covariant hot integrals are evaluated before the cold 
divergent integrals. Once the Lorentz invariance breaking hot loop energies are 
integrated out, the usual techniques of Feynman parametrization and dimen- 
sional regularization can be applied to perform the cold loop Lorentz invariant 
terms, as in vacuum. It is also worth-mentioning that all the hot corrections 
gives T^ dependance. The incorrect order of integration gives increasing order 
of T also (See Eq.(9)). At the two loop level, it can be seen from Eqs.(16) - (18) 
that the vertex type corrections to the virtual electrons vanish and the self en- 
ergy type corrections to the electron loops contribute at low temperatures. The 
dominant temperature corrections in this hot medium are the lowest order a, 
which ensures the convergence of perturbative series. Since the thermal effects 
are incorporated as perturbative effects, at sufficiently low energies {E << nie), 
the finite temperature corrections become negligible and we can always repro- 
duce the vacuum results by taking T=0 limit in our results. The presence of 
the statistical effects of photons modifies the vacuum polarization and hence the 
electron charge which leads to the changes (though small) in the electromag- 
netic properties of the hot medium even at low temperatures. The mass, the 
wavefunction, and the charge of electron are renormalized in the presence of the 
background heat bath and the calculated renormalization constants in QED are 
in the hot photon background up to the second order in a. These renormalized 
values can simply give the physical mass of electron and the effective charge and 
the wavefunction of the electron in such a background. Equations (25) and (26) 



give the estimation of the dielectric constant and the magnetic permeability of 
a hot medium. Both of these quantities deviate from unity whereas their values 
were unity at the one loop level. 

It may be noted that as p^ — > 0, Eq. (22) reduces to 



HtW^-^, (29) 

which corresponds to the dynamically generated thermal mass of photon in 
the hot background. This type of effect has been earlier observed for self- 
mass of gluon even at the one- loop level [6] . We further obtain the respective 
propagation vectors and frequencies when we take the limiting values for the 
longitudinal and transverse component of the vacuum polarization tensor in Eqs. 
(21) and (22). In particular, in Eq. (27), ^"^"^ represents the Debye screening 
length of electric force in such a medium. In this way, our results are not only 
an explicit proof of renormalizability of QED up to the two-loop level, they 
also estimate the temperature dependent modification in the electromagnetic 
properties of a medium. This helps to evaluate the decay rates and the scattering 
crossections [13] of particles in such media. This is a bit lengthy procedure, 
however it gives a good estimate of the background contribution. In our next 
projects, we would try estimate the renormalization constants of QED in hot 
and dense media at the two-loop level and would apply the same procedure in 
the presence of a constant magnetic field in the background. Such work is being 
done in other schemes [14] as well. 



6 APPENDIX 



The hot contribution at low temperature, comes from the photon background 
only. Eq. (12) for Fig. (2a) therefore simplifies to 

1 

(fc^ — TTi^ + ie)"^ [{p — k)^ — m? + ie\ [{k + V}^ — m? + ie\ 

where 



N^^, = 8 [(3to2 _ fc2 _ 2k.l)] k^{p - k), + k,{p - k)^ - g^,k.{p - k) + {^ - m^) {l^{p - k), + 
lu{p - k)f, - gf,J.{P - k)} + 2g^i,m^(m^ - k.l)], 

and 

10 



g-i^uip^ T)-^j ,Ji:s^::!)-L) / iii^iii-^« 



with 



|lM|l|nB(0 



2rp2 



TT^T 



Now integrating over the cold loop by using Feynman parametrization and 
dimensional regularization one gets Eq. (14). On the other hand, first do- 
ing Feynman parametrization, integrating by dimensional regularization, and 
simplifying one gets 

5-n^.(p, T) = ^ / dHS{P)nB{l) /o' dx /J= dy J^ dz 

J/ + ZJ+4m |H m2_2p.,y(y+^) ^^i;^ - 7 - in (^ -m^-2pjy(y+z) j /J' 

and 

c;-n-,(p,r) - ^/dA./|lM|lK(0/o rfx^f dy/„^dz{2(2y4+2yV+V;^_ 



-yz)|l|2 


(po-IpIa* 


ll+fr 


o+IpIm 
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+m2|l|(6y2+62/z-y+z) 


> 


-IpIm) 
— T — 


(po+IpIa*) 

T 

a_ 


|l|(24y2 + 24yz-18y-12z) 


(po-IpIa*) (po+IpIa*) 

2 ^ 

ai_ a_ 


+ Am^ 


1 + 1 




-6m2 " 


a_j_ a_ 


-12 


2(i-7-ln47r)-ln(a+)-ln(a_)]}, 


where 


a± = rr 


." ± %(9 


+ z)(j 


'o- + |P 


M 


. 











After a lengthy calculation, Eq. (38) finally simplifies to Eq. (15). 
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